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ABSTRACT 

We present a detailed photoionization model of the brightest knot of star formation 
in the blue compact dwarf galaxy Mrk 209. The model reproduces the intensities of 
main lines emitted by the ionized gas, resulting in a very good agreement between 
observed and predicted line temperatures and chemical abundances of the observed 
ionic species. 

The model has been calculated using the spectral energy distribution of a massive 
cluster of recent formation as the ionizing source. The features of Wolf-Rayet stars 
observed in the spectrum of the object, along with its ionizing properties, lead to dif- 
ferent solutions for the ages and characteristics of the ionizing stellar populations. The 
found solutions are compatible with either a composite population of two ionizing 
clusters or a continuous star formation. In the first case, a young cluster, with an age 
less than 3 Myr, would be responsible for most of the ionization properties while an 
older cluster, with either 3.2 or 4.2 Myr, would be responsible for the emission of the 
observed WR features in the spectrum of the knot. In the second case, the duration 
of the star formation episode is found to be 3.6 Myr. 

The addition of a certain amout of dust was required in order to reproduce cor- 
rectly the measured electron temperatures. This dust implies a gas/dust ratio of 
6.13xl0 -3 , identical to that in the Milky Way, and a fraction of absorbed photons 
of / = 0.58. The taking into account the dust grain physics combined with a thick- 
shell geometry solves the problem of the heating in this object and allows the fitting 
of the auroral lines, the line temperatures and the ionization structure of the nebula 
with a negligible presence of temperature fluctuations, in agreement with the most re- 
cent results found for this type of objects from the measurement of the Balmer jump 
temperature. 

An analysis of the ionization structure of the gas as well as the calculation of the 
main ionization correction factors for the unseen ionic stages is presented. 
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1 INTRODUCTION 

Blue compact dwarf galaxies (BCD) are objects whose spec- 
tra, as in the case of HII regions, are dominated by emis- 
sion lines from the gas ionized by very recent episodes of 
star formation (Sargent & Searle, 1970). Many intermediate 
redshift objects have properties similar to the BCD in the 
Local Universe and it is possible that they were common 
in the past and they have evolved to other kind of objects. 
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Therefore, the study of the star forming history in BCDs is 
of great importance. 

In general, BCDs are characterized by their compact 
aspect, very low metallicities, gas richness and blue colors 
(Kunth & Ostlin, 2000). These properties have made of them 
good candidates to host their first episodes of star formation. 
Nevertheless, recently, the detection in most of them of low 
surface-brightness elliptical haloes, or the presence of stars 
belonging to older populations , have caused this interpre- 
tation to be reconsidered. Nowadays, only few candidates 
remain controversial 

There are several works in the literature whose main 
aim is to study the weight of these older stellar populations 
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in BCDs and, hence, to find out the actual age and evo- 
lutionary status of these objects. Among these studies, in 
the sample of local objects, observations with enough spa- 
tial resolution to provide photometry of the individual stars 
have allowed, by means of colour-magnitude diagrams, to 
date some of the bursts (e.g. VIIZw403 by Schulte-Ladbeck 
ct al., 1998; IZwl8 by Aloisi et al., 2000). In all these works 
evidences for stellar populations older than 500 Myr have 
been found. Besides, old stellar low-surface brightness com- 
ponents have been detected in these local objects by study- 
ing the radial light distribution in the optical (e.g. Gil de 
Paz et al., 2003) and the near infrared (e.g. Noeske et al., 
2005). 

For more distant galaxies, the study of the stellar popu- 
lation is limited to the interpretation of the integrated spec- 
trophotometry properties by means of empirical and evo- 
lutionary population synthesis techniques. The results from 
this kind of works (e.g. Raimann et al., 2000) show that 
most blue compact dwarves are age-composite stellar sys- 
tems, with discontinuous star-forming histories. 

Nevertheless, there is still a lack of spatially resolved 
studies of the properties of the ionizing stellar populations 
and the ionized surrounding gas in the different knots of 
star formation. This work is usually approached by means 
of photoionization models, trying to reproduce the measured 
fluxes and equivalent widths of the lines emitted by the ion- 
ized gas, under different assumptions about its geometry, 
metal content and ionizing incident radiation. Several im- 
provements have been produced that increase our confidence 
in the results from this kind of models: inclusion of blan- 
keting effects in O and WR expanding stellar atmosphere 
models (Pauldrach et al., 2001; Smith et al., 2002) and re- 
cent computation of NLTE blanketed plane-paralel models 
(Lanz & Hubeny, 2003). Yet, model computed evolutionary 
sequences show important disagreements with observations. 
Among others, the [Om]/H/3 vs [Oll]/H/3 and [Om]/H/3 vs 
[OI]/H/3 relations are not well reproduced by evolutionary 
model sequences in the sense that predicted collisionally ex- 
cited lines result too weak compared to observations, which 
would require additional heating mechanisms or the pres- 
ence of chemical inhomogeneities (Stasiriska & Izotov, 2003) 

These facts seem to point out to a still incomplete de- 
scription of the ionization structure of the emitting nebulae 
which would improve with the determination of more than 
one line electron temperature. This has been attempted re- 
cently by several authors (e.g. Garnett, 1989; Perez-Montero 
& Diaz, 2003, hereinafter PMD03; Hagele et al. 2006) by en- 
larging the spectral sampling to include the lines of [On] at 
AA 7319,30 A and [Sin] at AA 9069, 9532 A. The far-red [Sm] 
emission lines can provide better estimates of the ionization 
parameter of the gas, via the [Slll]/[Sn] quotient (Diaz et al., 
1991), and the equivalent effective temperature of the ion- 
izing clusters, via the rj parameter (Vflchez & Pagel, 1988), 
than the strong widely used optical lines of oxygen, mainly 
due to their lower dependence on other quantities, such as 
age or metallicity. These parameters are key factors to pro- 
vide adequate input constraints for photoionization models. 

Although the detection and analysis of the near-infrared 
emission lines is hampered by the strong and variable ab- 
sorption features of the Earth's atmosphere, a thorough at- 
mospheric extinction correction and a good spectral reso- 



lution aid to remove the sky background providing a good 
level of confidence in the measurement of the far red [Sin] 
lines. The number of Hn galaxies observed both in the op- 
tical and the far-red spectral range, up to 1/i, is still scarce 
but it is steadily increasing. 

In this work, we apply a detailed photoionization model 
to the brightest knot of the BCD Mrk209, for which both 
wide range spectroscopic observations and photometric data 
exist. 

In the next section we describe the properties of Mrk209 
as given in the literature, focusing on the main burst of star 
formation, measuring their photometric properties and the 
gas physical conditions. In section 3 we describe the pho- 
toionization model we have used with the inclusion of the 
input stellar ionizing population and the effects of dust. Sec- 
tion 4 presents and discusses the model results in comparison 
with observations. Finally, the conclusions of this work are 
presented in section 5. 



2 PROPERTIES OF THE OBJECT 

Mrk 209, also known as IZw36 and UGCA 281, is a blue com- 
pact dwarf galaxy that has been extensively studied. The 
work by Viallefond & Thuan (1983) revealed a core-halo 
structure of the diffuse neutral gas, with the star-forming 
regions located near the core, but slightly shifted with re- 
spect to the peak in Hi column density. They measured its 
virial mass which is about 6 times that in Hi, which could 
imply a large contribution of dark matter. Loose & Thuan 
(1986) discovered an extended, elliptical background sheet 
underlying the compact, actively star forming core and they 
classified this galaxy as iE, corresponding to blue compact 
dwarves with an irregular burst of star formation in its core, 
surrounded by regular elliptical isophotes. 

Although Fanelli et al. (1988), from an IUE UV spectral 
study, pointed to Mrk209 being a young galaxy undergoing 
its first episode of star formation, later works have ruled 
out this conjecture. The red colours observed by Deharveng 
(1994) and by Papaderos et al. (1996) require the presence of 
old stellar populations in the elliptical host galaxy. Schulte- 
Ladbeck et al. (2001) were able, by means of HST observa- 
tions with FOC in the optical and NICMOS in the near-IR, 
to build a colour-magnitude diagram identifying individual 
stars of the red giant branch (RGB), some of them with ages 
between 1 and 2 Gyr. These authors claim for an almost- 
continuous star formation history for this galaxy with inac- 
tive periods never longer than 100 Myr. 

Gil de Paz et al. (2003; hereinafter GMP03) have pro- 
duced photometric observations in Ha and the B and R 
bands. They find that the Ha emission peaks at a region 
towards the west of the galaxy, accompanied by two weaker 
star forming regions towards the east. The weakest of these 
regions lacks near infrared emission (Noeske et al., 2005) 
implying that the dominant population is still very young. 

Regarding the brightest knot of star formation, Dehar- 
veng et al. (1994) showed that their young stars have an 
age less than 12 Myr. Later evolutionary synthesis models 
by Mas-Hesse & Kunth (1999) have lowered this upper limit 
further to only w 2.7 Myr. Besides, Izotov et al. (1997; here- 
inafter ITL97) and Schaerer et al. (1999) discovered WR 
features in this knot. 
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Table 1. Properties of the main burst of star formation in Mrk 
209. 



Distance (Mpc) 


6.6±1.1 


F(Ha) a 


14.97±0.01 


Z b 


0.0036 




0.06 


B (mag) c 


16.58±0.13 


B - R (mag) 


0.28±0.18 


EW(H/3) (A) 


206 


EW(4650 A) (A) 
log L(4650 A) d 


5.01±0.10 


37.35±0.13 


EW(5808 A) (A) 


1.74±0.15 


log L(5808 A) e 


36.62±0.16 



a in units of 10 13 erg ■ s 1 -cm 2 

6 in mass, scaled to the oxygen solar abundance from Allende- 
Pricto, 2001 

c not corrected for reddening 
d in units of erg ■ s _1 



Table 2. Emission line intensities, relative to H/3= 100, with their 
corresponding errors, observed and the predicted by our model 
and the model in PMD03. 



Emission line 


Observed 


This model 


PMD03 


3727 A [Oil]" 


71.9±0.2 


69.4 


68.7 


3868 A [Nem] a 


45.7±0.1 


46.0 




4072 A [Sn] a 


1.2±0.1 


1.2 


_ 


4363 A [Om] a 


12.7±0.1 


12.4 


9.6 


4658 A [Fem] a 


0.3±0.1 


0.4 




4959 A [Om] a 


196.0±0.3 


190.1 


200.8 


5007 A [Om] a 


554.3±0.8 


574.2 


604.4 


6584 A [Nn] b 


3.6±0.1 


3.6 




6725 A [Sn] b 


10.6±0.4 


10.2 


6.8 


6312 A [Sin] 6 


1.9±0.1 


1.9 


1.1 


7137 A [Ann] 6 


5.5±0.1 


5.6 




7325 A [On] 6 


1.7±0.2 


2.6 


2.4 


9069 A [Sm] 6 


12.2±1.2 


12.1 


8.7 


10.5 j± [Siv] c 


42.6±7.6 


16.5 




18.7 n [Sm] c 


17.4±4.5 


9.7 





Spcctrophotometric observations in the range between 
3000 and 7800 A from ITL97 and between 7000 and 9800 A 
(PMD03) exist for this burst. The slit width used in both 
works is identical and the agreement of the relative fluxes 
in the coincident spectral range is good; therefore, we can 
be confident about the spatial coincidence of both sets of 
observations. The intensities and equivalent widths of the 
observed WR features were measured by Guseva et al. (2000; 
hereinafter GOO). 

We have collected the main photometric and spectro- 
scopic properties of the brightest knot of star formation in 
Table 1. Regarding the spectroscopic data, WR feature in- 
tensities and equivalent widths (GOO) refer to a slit width 
of 2 arcsec, while the emission line fluxes, H/3 equivalent 
width and reddening constant c(H/3) (ITL97; PMD03) re- 
fer to a slit width of 1.5 arcsec. Besides, we have compiled 
ISO observations of this object (Nollenberg et al., 2002) that 
account for line intensities of 4.05 [s, Bra, 10.5 \jl [Siv] and 
18.7 fj, [Sill], in the mid infrared. Nevertheless the apertures 
used for these observations are much larger that those used 
in the optical and the near infrared and, althought it is not 
expected to affect the fluxes of the lines of the ions of high 
excitation it could do it to the hydrogen recombination line. 

Luminosities have been derived assuming a distance of 
6.6±1.1 Mpc (Schulte-Ladbeck et al. 2001). The quoted er- 
rors for luminosities include the uncertainty in the distance 
to the galaxy. Regarding the photometric data, we have se- 
lected the area of the Ha image from GMP03 containing 
the same flux as in GOO corrected for aperture, and we have 
re-measured the B and R magnitudes in the corresponding 
images. The selected region has a mean radius of 72±18 pc. 
The reddening corrected emission line intensities, relative to 
H/3 = 100, are given in Table 2. 

Using these emission line data we have recalculated the 
electron temperatures and density, and the ionic chemical 
abundances. We have used the same algorithm described 
in PMD03: the five-level statistical equilibrium model in 
the TEMDEN and IONIC tasks of the NEBULAR pack- 
age in STSDAS, based on the FIVEL program developed 



a from Izotov et al., 1997 

b from Perez-Montero & Dfaz, 2003 

c from Nollenberg et al., 2002 

by De Robertis et al. (1987) and improved by Shaw & 
Dufour (1995). We have used the transition probabilities 
and collision strengths included in the photoionization code 
CLOUDY (version 96.0; Ferland et al., 1998), except in the 
case of + for which we have used the transition probabil- 
ities from Zeippen (1982) and the collision strengths from 
Pradham (1976), which offer more reliable results for nebu- 
lar diagnostics (Wang et al., 2004) for this species. 



3 MODEL DESCRIPTION 

The model here described simulates the properties of the 
ionized gas and the stellar ionizing population of the bright- 
est knot of star formation in Mrk 209. It has been calculated 
using the photoionization code CLOUDY (Version 96.0, Fer- 
land et al., 1998). This model is characterized by a set of in- 
put parameters including the ionizing continuum, the nebu- 
lar geometry, the gas density and chemical abundances. The 
main properties of the best fitted model are summarized in 
Table 3. 

3.1 The ionizing stellar population 

The WR phase in the evolution of massive stars is a short 
and transient stage, and therefore the presence of WR fea- 
tures in the spectrum of a star forming region allows a dating 
of the ionizing stellar population. The emission features usu- 
ally observed consist of two broad bumps. The first, usually 
called as blue is at 4650 A and is a blend of the Nv A 4605, 
4620 A, Nm A 4634, 4640 A, Cm/Civ A 4650, 4658 A and 
Hen A 4686 A broad WR lines. These are emitted mainly 
by WN and early WC stars. Besides, early WN stars emit 
broad emission features of Civ at A 5808 A, known as the 
red bump. Although the number of WR stars is expected to 
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Table 3. Properties of the modelled Hll region. 



3.3 Dust content 



Observed This model 



EW(H/3) 206 A 210 A 

log Q(H) a 51.75±0.13 51.96 6 

log U -2.3 ± 0.1 -2.22 

filling factor 0.077 

Inner radius 30 pc 

Stromgren radius 72±18 pc 48 pc 

Absorption factor / 0.58 

Dust-to-gas ratio 6.13 • 10~ 3 

C(H/3) 0.06 0.19 



a in photons ■ s 1 

b not absorbed by dust 



be much lower in relation to the number of O stars at low 
metallicites, they have been detected in the spectra of some 
blue compact dwarfs, including Mrk209 (ITL97, Schaerer et 
al., 1997) implying the presence of stars in the burst with an 
age between 1 and 5 Myr, that is the age of the WR stars. 

We have calculated the spectral energy distributions of 
ionizing star clusters in this range of ages using evolutionary 
synthesis techniques (Starburst 99; Leitherer et al., 1999) 
with the inclusion of the implementation by Smith et al. 
(2002) using the model atmospheres from Pauldrach et al. 
(2001) and Hillier & Miller (1998) for O and WR stars re- 
spectively which include the effects of blanketing. The used 
initial mass function is a Salpeter's one, with masses between 
0.8 and 120 Mq. The supernova cut-off mass has been set to 
8 Mq and the black-hole cut-off mass is at 120 Mq. We have 
considered a metallicity for the models equal to 0.004, 1/20 
times the solar value, which is very close to the metallicity 
derived for this object, equal to 0.0036. We have considered 
models with both standard and high mass loss to explore 
which models fit better the observed properties of the ob- 
ject. 



3.2 Nebular geometry 

One of the input parameters of the models is the inner radius 
of the ionized region. In the most successful model, which is 
the one presented here, the ionization front is located at a 
distance of 30 pc of the ionizing source. Since the thickness 
of the Hll region reaches approximately 18 pc, the resulting 
nebular geometry is a thick shell. The total radius, therefore, 
is not far from that measured on the Ha image around the 
brightest knot, 72±18 pc, encompassing the same flux mea- 
sured in the slit and corrected for aperture effects by GOO. 
This geometry, along with a filling factor of 0.077 leads to 
a value of the spherical ionization parameter log U = -2.22, 
coincident with the value -2.3±0.1 calculated by PMD03. 
This model considers a constant density of 190 particles per 
cm 3 , which is the value obtained from the quotient of fluxes 
of the lines of [Sn] . 



We have added a certain amount of dust in order to fit 
correctly the measured electron temperatures. The pres- 
ence of solid grains within the ionized gas have some con- 
sequences on the physical conditions of the nebula that 
should not be neglected, including the depletion of metals 
in gaseous form onto grains or the absorption of energetic 
radiation. The heating of the dust can affect, hence, the 
equilibrium between the cooling and heating of the gas and, 
thus, the electron temperature inside the nebula. The com- 
mand PGRAINS, included in the last version of Cloudy, 
implements the presence of dust grains as described in van 
Hoof et al. (2001). This command resolves the grain size dis- 
tribution into several bins and treats each bin with its cor- 
responding temperature, potential and drift velocity, which 
depend on the diameter of each grain. We have assumed the 
default grain properties of Cloudy 96, which has, essentially, 
the properties of the interstellar medium and follows a MRN 
(Mathis, Rumpl & Norsieck 1997) grains size distribution. 



4 RESULTS AND DISCUSSION 

We find that the model that better fits the observed line 
intensities, the derived ionic abundances and the other ob- 
servational properties of this knot is well represented by an 
instantaneous burst of 1.0 Myr. The inclusion of another 
burst in the WR phase in order to explain the observed fea- 
tures does not affect the ionization structure of the nebula 
and the obtained results. This result is equivalent assuming 
a constant star formation rate for the ionizing burst in such 
a way that its present age was that in the WR phase. The 
predicted emission line intensities relative to H/3 = 100 are 
presented in column 3 of Table 2 and can be compared with 
the observed ones listed in column 2, and with those pre- 
dicted by the best-fitting model in PMD03. 

We can see all the relevant emission lines are bet- 
ter fitted by the present model than by that presented in 
PMD03 in which the ionizing radiation was represented by 
the CoStar atmosphere of a single massive star and a plane- 
parallel geometry was asumed. Using that model, it was not 
possible to reproduce the observed T([Om]) due to the un- 
derestimate of the auroral line at 4363 A. On the other hand, 
an incorrect simultaneous fitting of the strong lines of [On] , 
[Om], [Sn] and [Sin] caused an overestimate of the r}' pa- 
rameter, defined as 

, _ [On]3727/[Om](4959 + 5007) 
71 ~ ([S'ii](6717 + 6731)/[S , m](9069 + 9532) 

(Vilchez & Pagel, 1988), which is indicative of the temper- 
ature of the radiation field. In the present model, by means 
of including a more reallistic ionizing spectral energy distri- 
bution, derived from the evolutionary synthesis of the ion- 
izing cluster, dust physics and an appropriate geometry, we 
have been able to reproduce adequately the main emission 
line intensities. Nevertheless, the model is not able to repro- 
duce the mid-infrared lines observed with ISO, obtaining 
values for both [Sin] 18.7/x and [Siv] lower than observed by 
60% and 40% respectively. This illustrates the problem that 
photoionization models have in general to correctly repro- 
duce the intensities of sulphur emission lines in this spectral 
range. In Figure 1 we represent the quotients of the lines 
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# SpitzerBCDs (Wu etal. 2006) 

□ ISO on Mrk209 (Nollenberg et al. 2002) 

X This model for Mrk209 
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log U = -2.0 
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IIZw40 
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Spitzer BCDs (Wu et al. 2006) 

□ ISO on Mrk209 (Nollenberg et al. 2002) 

X This model for Mrk209 




- UM46li 

4 

3 T 



IIZW40 TIZw18 
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log ([Oil] 3727/[OIII] 4959+5007) 



Figure 1. Diagnostic diagrams of the quotients of [Siv] 10.51 fi / [Sin] (9069+9532 A), at left and of [Sill] 18.67 / [Sll] (6717+6731 
A), at right, versus the quotient of [On] 3727 A/ [Om] (4959+5007 A), for a sample of Blue Compact Dwarf Galaxies observed both 
in the near and the mid infrared and for a sequence of models covering the physical condicitons of ionized gaseous nebulae. All model 
sequences, including the model presented here, fail to predict correctly the intensities of the sulphur emission lines in the mid-IR. 




2 4 60 2 4 6 

Age (Myr) Age (Myr) 

Figure 2. The colour index B-R predicted by evolutionary synthesis models for both standard and high mass loss stellar evolution 
and for instantaneous (thick solid line) and continuous (dashed line) star formation rates. The bands represent the values measured in 
the brightest star forming region in the photometric images. 



of [Siv] 10.51 m / [Sm] (9069+9532 A), at left and of [Sm] 
18.67 fi I [Sn] (6717+6731 A), at right, versus the quotient 
of [On] 3727 A/ [Om] (4959+5007 A), for a sample of BCDs 
observed both in the mid-IR with Spitzer (IZwl8, UM461 
and IIZw40; Wu et al., 2006) and ISO (Mrk209; Nollenberg 
et al., 2002) and the near IR: UM461, IIZw40 and Mrk209 
(PMD03), IZwl8 (Garnett, 1992). The comparison with the 
model presented here and a sequence of CLOUDY models 
simulating a very young cluster ionizing nebulae of different 
physical conditions show serious discrepancies with observa- 
tions for the prediction of the fluxes of the sulphur lines at 
mid-IR. 



4.1 Inner extinction 

The inclusion of grain physics helps to reproduce the in- 
tensity of the [Om] at 4363 A, which depends exponen- 



tially on electron temperature. This shows the relevance of 
the presence of dust in the nebula. The required amount of 
dust implies a dust-to-mass ratio equal to 6.13T0 -3 , almost 
identical to that found for the Milky Way (Spitzer 1978). 
This amount of dust for the column density predicted by 
the model, implies a visual extinction of Av = 0.39 mag 
and would correspond to a logarithmic extinction at H/3 of 
0.19. This is higher than derived from the Balmer decrement 
(C(H/3 = 0.00; ITL97, C(H/3 = 0.06; GOO), but consistent 
with what is expected when comparing the computed ion- 
izing star cluster colour index (B-R) with the observed one. 
Assuming that the presence of any older stellar population 
has not a relevant weight for this emitting knot, we can cal- 
culate the effect of the extinction by the dust present in the 
nebula on the synthesized (B-R) colour. The predicted (B- 
R) colour index for the synthesized clusters as function of 
age is shown in Figure 2. The average colour predicted by 
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Z = 0.004, standard mass loss 
WR at 4650 A WR at 5808 A 



Z = 0.004, high mass loss 
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Figure 3. Comparison between the intensities and equivalent widths of WR red and blue bumps observed in Mrk209 and those 
predicted by the models. At left, models with standard mass loss and at right, with high mass loss. Both have been calculated for a 
metallicity Z = 0.004. The clear bands stand for the observed values with their errors as published in GOO. The dark bands stand for 
the observed values once corrected for nebular continuum, in the case of EW and the fraction of absorbed ionizing photons, in the case 
of relative intensities. Models for an instantaneous burst are plotted as a solid thin black line. Models for a continous star formation 
are plotted with a dashed line. Finally, the solid thick line represents models of a composite population with an instantaneous burst 
with an age of 1 Myr, resulting in the expected EW(H/3). 
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Figure 4. Prediction of the models for the equivalent width of H/3 for both standard and high mass loss and for instantaneous (thick 
solid line) and continuous (dashed line) star formation rate. The horizontal line represents the observed EW(H/3, once substracted the 
nebular continuum and corrected for the amount of absorbed ionizing photons as predicted by our model. 



instantaneous bursts of star formation for the ages between 
2 and 5 Myr assumed in our models is around -0.05. This 
is bluer than the colour index measured for the knot which 
is B-R = 0.28±0.18. However, according to our model, the 
predicted colour should be reddened by 0.24 magnitudes for 
the required amount of dust, which would make it consistent 
with observations. 

The effect of dust has to be considered when calculating 
the number of ionizing photons reaching the gas. In fact, 



only a fraction / of the number of ionizing photons emitted 
by the cluster will ionize the gas 



Q(H) = f-Q (H) 

For the amount of dust required to reproduce the emission 
line intensities in our model, we obtain / ~ 0.58. 
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Figure 5. Comparison between the intensities and equivalent widths of the WR red and blue bumps observed and aperture corrected 
in Mrk209 and those predicted by the models. At left, models with standard mass loss and at right, with high mass loss. Both models 
have been calculated for a metallicity Z=0.004. The bands stand for the same quantities represented in Figure 3, but now taking into 
account the aperture correction. Models for an instantaneous burst are plotted as a solid thin black line. The solid thick line represent 
models of a composite population with an instantaneous burst with an age of 1 Myr, resulting in the expected EW(H/3) in the slit 
(443A) and that assuming the aperture correction only for the 1.0 Myr population (514A) 



4.2 ionizing stellar populations 

The comparison between the observed intensities and equiv- 
alent widths of the blue and red bumps and those predicted 
by the models is shown in Figure 3. The light- coloured bands 
represent the intensities and the equivalent widths of the 
WR features with their corresponding errors as measured by 
GOO inside the slit. The dark-coloured bands represent the 
same quantities corrected for different effects: In the case of 
equivalent widths, for the contribution of the nebular contin- 
uum that is 20% of the ionizing continuum at the wavelenght 
of the blue bump, and 25% at that of the red bump as pre- 
dicted by our model. For the relative fluxes, the correction 
is due to the absorption of ionizing photons by dust, affect- 
ing the intensity of H/3. Regarding models the solid thin line 
represents the predictions by a model with a unique instan- 
taneous burst whose age is that shown by the X axes, and 
the dashed line the prediction for a burst with continuous 
star formation. Agreement between the observed and pre- 
dicted relative intensities and equivalent widths of the WR 
features is seen only for the models with a high mass loss 
rate and for an instantaneous burst of aproximately 3.0 Myr. 

Another important constraint to the model is the equiv- 
alent width of H/3, which is an estimator of the age of the 
ionizing burst. We consider as well negligible the contribu- 
tion of the underlying stellar population in this knot with 
a very high value of EW(H/3). Again the value of 206 A re- 
ported by GOO has to be corrected for the contribution of 
nebular continuum and absorption of ionizing photons by 
dust. These effects increase the value of EW(H/3) produced 
by the cluster to 443 A. This value is much higher than pre- 



dicted by Starburst 99 models for a cluster of 3.0 Myr and 
high mass loss rate as it can be seen in Figure 4. However, 
this disagreement would be solved if a younger population 
is present. We have therefore explored this possibility by 
assuming that the ionizing population is composed by two 
different bursts, one responsible for the WR emission and 
another that we have taken of 1.0 Myr that would provide 
the necessary excess of ionizing photons. The predicted val- 
ues of WR equivalent widths and fluxes for this composite 
population whose EW(H/3) is that observed and corrected, 
are shown by solid thick lines in Figure 3. It can be seen 
that no agreement is now found for the blue bump although 
the red bump values are well reproduced. 

Part of this disagreement could be due to aperture ef- 
fects. GOO point to a factor of 3.1 between the H/3 flux mea- 
sured inside the slit and the total H/3 emission. If we assume 
that correction and consider the WR stars to be concen- 
trated inside the slit, we obtain the fluxes and EW of the 
WR features shown in Figure 5. The measured fluxes and 
equivalent widths have been divided by the aperture correc- 
tion factor, as a consequence of the increment of the emission 
of H/3 and its continuum and assumming it affects similarly 
to both quantities. Again, we plot the same quantities cor- 
rected for the presence of the nebular continuum and the 
absorption of ionizing photons by dust. The predictions of 
the models are plotted with a solid thin line for a unique 
burst of star formation and a dashed line for continuous 
star formation. The solid thick line corresponds to a com- 
posite population with an instantaneous burst of the age 
showed at the X axis, plus a population of 1.0 Myr, the two 
of them together providing the expected EW(H/3). Now, we 
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find agreement for both blue and red bump features in the 
high mass loss rate case assuming either a constant star for- 
mation rate at an age of 3.6 Myr or a populations composed 
of a young cluster of 1.0 Myr and an older one with age 3.2 
or 4.2 Myr which contains the WR stars. In the case of the 
continuous star formation rate, the agreement occurs for an 
age slightly higher than that compatible for the expected 
EW(H/3) , but probably within the errors. 

We can calculate the mass of ionizing clusters in each 
case, using the number of ionizing photons, which can be 
obtained from the flux of Ha, from the relation: 

logQ(-ff) = 11.86 + log L(Ha) 

Taking into account the fraction, f, of photons absorbed 
by dust, the resulting number of ionizing photons is 



10 



r,1.99±0.13 



Neglecting the contribution to the number of 



ionizing photons by the underlying population : 

Q{H) = M yc ■ Q(H) yc + M WR ■ Q(H) WR 

where M yc and Mwr are the masses of the clusters of 1.0 
Myr and in the WR phase, respectively, and Q(Ha) yc and 
Q(Ha)wR, their corresponding numbers of hydrogen ion- 
izing photons per unit mass. The relative contributions of 
each of the clusters to the total ionizing mass can be derived 
from the observed equivalent width of H/3 using the equation 
below: 

EW(Hf3) obs - EW{Hf3) yc = L c (Hf3) WR = 
EW(Hf3) WR - EW(Hf3) ob3 L c (Hf3) yc 

_ L c (Hf3)wR/m M WR 
~ Lc(Hf3) yc /m M yc 

where EW(H/3) 6 S stands for the observed value of the 
equivalent width of H/3, EW(H/3)„ C , EW(R/3)wr are the 
predicted equivalent widths for each of the clusters and 
Lc(H/3) !/c /m, Lc(H/3)ivij/m are the model predicted con- 
tinuum luminosities per unit mass at the H/3 wavelength for 
each of the clusters. 

The characteristics of the model ionizing clusters are 
listed in Table 4. The luminosities for the single ionizing 
clusters are given per unit solar mass. Using these values and 
the expresions above, it is possible to calculate the cluster 
masses for each assumed case. For the composite population 
with clusters 1 and 3.2 Myr old, masses of 8.7 x 10 4 and 2.4 
x 10 4 M Q are found. For the 1 and 4.2 Myr old clusters, the 
corresponding values are 9.3 x 10 4 and 3.3 x 10 4 Mq. This 
gives a total mass for the ionizing population of 1.11 x 10 5 
and 1.26 x 10 5 M respectively. In both cases the younger 
cluster makes up more than 70 % of this total mass and 
dominates the ionization providing more than 80 % of the 
hydrogen ionizing photons and about twice the continuum 
luminosity at H/3 than the cluster containing most of the 
WR stars. The observed luminosities of the blue and red 
bumps are better reproduced by the combination with the 
younger WR cluster, as well as the equivalent widths of the 
WR features. On the other hand, a larger cluster mass value 
is found for the continuous star formation scenario: 6.0 x 
10 5 Mq with the masses derived from the WR blue and red 
luminosities being barely consistent with observations. 

The found ranges for the mass of the starburst lie in a 
regime where stochastic fluctuations of the number of WR 
stars in relation to the expected value can still be important. 
In fact, the dispersion calculated by Monte Carlo simulations 



(Cervino et al., 2002) for the total mass of the clusters reach 
0.01 in F(WR)/F(H/3) and 1 A in the EW(WR) for both the 
blue and red bumps. This can affect the expected age of the 
cluster in the WR phase, but the condition of fitting at same 
time the equivalent width of H/3 constrains the possible age 
to only 1 Myr around the obtained value. Nevertheless, wc 
do not find serious discrepancies between the prediction of 
both scenarios, a composite population and a continous star 
formation about the number of WR stars, which is approxi- 
mately 4 or 5 WC stars and 1 WN star, and those obtained 
from comparing the luminosities found for both the blue and 
the red bump, with the mean luminosities of WC and WN 
stars. 



4.3 Line temperatures, ionic abundances and 
ionization correction factors 

The best-fitted model for the emission lines have been ob- 
tained by varying the properties of the geometry and the 
metal content of the gas for the obtained stellar energy dis- 
tribution. Here we summarize for each ionic species the val- 
ues of electron temperature and ionic abundances predicted 
by this model. In figure 6 we represent the radial profiles of 
the abundances of each ionic species in the best-fitted model 
along with the electron temperatures in order to illustrate 
the most representative temperature in the zones where each 
ionic species stays. 

In order to calculate averages of the electron temper- 
ature, ionic fractions and ionization correction factors, wc 
have considered volumetric integrations which fit better the 
observed values for models with a spherical geometry (Lurid- 
iana et al., 2002) 



Rs 



N(X l )N e dR 




N{Xi)N e dR 



The values obtained by the best-fitted model, along 
with the values obtained from the observations are tabu- 
lated in Table 5. 

Oxygen. Both [On] and [Oni] electron temperatures are 
relatively well reproduced by the model. While T([Oni]) re- 
sults 2% lower than the measured value, T([0n]) is overes- 
timated by about 10%. However, the lines involved in the 
determination of this temperature are widely separated and 
the measured error, mostly due to uncertainties in the red- 
dening correction, is large. Using the measured [Oni] tem- 
perature and the grid of models from PMD03, assuming the 
electron density obtained from [Sn] lines, we would obtain 
a value for T([Oll]) equal to 14300 K leading to a value 
for 12+log(0+/H + ) equal to 6.89±0.11, which is more co- 
incident with the result predicted by the model. Although 
the amount of + is not especially relevant for the total 
abundance of oxygen for this object, due to its high degree 
of excitation, this amount becomes more important for the 
calculation of the ionization correction factors of some other 
species through the + /0 quotient. The model gives as well 
a value of T([Oni]) via the recombination lines allowing to 
estimate the fluctuations of temperature inside the nebula, 
which could occasionate an underestimate of the chemical 
abundance (Peimbert, 1967). The value of 16100 K for this 
temperature leads to a t 2 = 0.004, which is consistent with 
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Table 4. Characteristics of the ionizing clusters as derived using the calculated number of ionizing photons and equivalent width of H/3. 





Young 


WR 


Composite 


("Ion tin n mis 


O VmprvpH 




elimtor 

KjI Liu Lj^_.± 


cluster 


Population 


SF 


( cnrrpctpH i 




(per M Q ) 


(per M ) 






(per M Q ) 




Age (Myr) 


1 


3.2 


4.2 


1, 3.2 


1, 4.2 


3.6 




logL(H/3) a 


34.62 


34.47 


34.16 


39.64 


39.61 


33.85 


39.63 ± 0.13 


logL(WR, 4650) a 


31.27 


33.07 


32.70 


37.46 


37.22 


31.60 


37.57 ± 0.13 


logL(WR, 5808)" 




32.71 


31.81 


36.57 


36.32 


31.11 


36.74 ± 0.16 


EW(H/3) (A) 


555 


220 


167 


443 


443 


427 


443 


EW(4650) (A) 


0.22 


2.19 


4.84 


2.40 


1.40 


2.03 


1.94 ± 0.04 


EW(5808) (A) 




7.63 


1.19 


0.72 


0.13 


1.23 


0.70 ± 0.06 



a in units of erg • s 



Table 5. Observed and model-predicted electron temperatures 
and ionic abundances. 





Measured 


Model-predicted 


T([On]) (K) 
12+log(0 + /H+) 

T([Om]) (K) 
12+log(0 2+ /#+) 

12+log(0/H) 


12400±1100 
7.10±0.11 
16200±100 
7.68±0.01 


13930 
6.88 

15850 
7.74 
7.80 


T([Nn]) (K) 
12+log(iV+/H+) 
ICF(N+) 
12+log(N/H) 


5.64±0.11 


13520 
5.50 
9.38 
6.47 


T([Neni]) (K) 
12+rog(Ae 2+ /#+) 
ICF(Ne 2+ ) 
12+log(Ne/H) 


7.01±0.01 


15600 
6.98 
1.05 
7.00 


T([Sn]) (K) 
12+log(S'+/H+) 

T([Sm]) (K) 
12+log(5 ,2+ /H+) 
12+log(S 3 +/# + ) 

12+log(S 3+ /# + )° 
ICF(S++S 2 +) 
12+log(S/H) 


13300±1600 
5.12±0.11 

15900±1700 
5.86±0.10 
5.93±0.20 
5.68±0.23 


13270 
5.05 

15600 
5.94 
5.52 
5.52 
1.31 
6.11 


T([Arni]) (K) 
12+log(Ar 2 +/ii'+) 
ICF(Ar 2+ ) 

T([Anv]) (K) 
12+\og(Ar 3+ /H+) 
ICF(Ar 2+ +Ar 3 +) 

12+log(Ar/H) 


5.31±0.09 
4.49±0.05 


15350 
5.32 
1.22 

15995 
4.59 
1.02 
5.40 


T([Fem]) (K) 

12+log(Fe 2+ /# + ) 
ICF(Fe 2+ ) 
12+log(Fe/H) 


4.92±0.22 


14450 
5.01 
7.59 
5.89 


Assuming a constant S 3 +/S 2 + ratio. See text for details. 



the measurements for Blue Compact Dwarf Galaxies made 
using the temperature of the Balmer jump (Guseva et al., 
2006: Hagele et al., 2006), which point to values of t 2 negli- 
gible in these kind of objects. 

The total oxygen abundance is calculated under the as- 
sumption that the fractions of neutral hydrogen and oxygen 
are equal: 

O ~ H 

due to the charge-exchange reaction + + H — ► O + H + 
(see Osterbrock, 1989), and this allows to calculate O/H as: 

o _ o+ + o 2 + 

H~ H+ 

In our model the neutral fraction of both elements are al- 
most coincident: 0.021 for H°/H and 0.019 for 0°/0. 

Nitrogen. T([Nll]) can be determined directly when the au- 
roral line, at 5755 A is measured. Nevertheless in the most 
objects, as in the case of Mrk209, this line is not detected 
with enough signal-to-noise ratio. Therefore, the ionic abun- 
dance of N + is calculated under the assumption that T([Nn]) 
« T([On]), since both ions lie in the low excitaton zone of 
the nebula. In our model, the electron temperature associ- 
ated to the zone of N + is slightly lower than that of + and, 
in fact, is closer to the value of T([Sll]). The N+/H+ ionic 
abundance calculated with T([0n]) as derived from the mea- 
sured T([Om]) using the grid by PMD03 is 12+log(N+/H+) 
= 5.51±0.02, almost identical to the value predicted by our 
model. 

The total abundance of nitrogen is calculated under the 
assumption that 

N+ _ Q+ 

~W ~ IT 

For our model, this approximation yields a value of ICF(N + ) 
equal to 8.24, confirming the results predicted by the mod- 
els in PMD03. The predicted value is slightly higher than 
the value obtained for the hot model from Mathis & Rosa 
(1991), which predicts ICF(N+) = 6.44. 

Neon. The ionic abundance of Ne 2+ is calculated assuming 
that T([Nem]) w T([Oni]), what seems well justified in view 
of the results of our model. The total abundance of neon is 
calculated assuming that 



© 2006 RAS, MNRAS 000, 1-?? 



10 E. Perez- Montero & A.I. Diaz 




0.8 



33 36 



39 42 
radius (pc) 





1 


i 1 






T[SIII] T 




S + - 


- 




T[SII]\A/ 


/ ~ 








- 


- 


s 2+ 




\ - 




















j , 1 , 1 1 r 


1 \. 





1 1 ' 


1 1 


i 1 










o . 






Ai 




A 2+ 








Ar 






A 3+ 








^--~^Ar 












Ar + , 




i 


i , — i — 







33 36 



39 42 
radius (pc) 



6000 




39 42 
radius (pc) 



Figure 6. Radial profiles of the relative abundances of the most representative ionic species: oxygen, sulphur, nitrogen, neon, argon 
and iron. The radius indicates the distance to the ionizing cluster, with the ionization front located at 30 pc.The ionic fractions are 
marked on the left y axes. The radial profile of the electron temperature is plotted for all the species along with their available measured 
values, which have been located at a distance that is a weighted average of the ionic fraction of the corresponding element. Its scale is 
marked on the right Y axes in all the plots. 



Ne 2+ _ Q 2+ 
Ne ~ ~0~ 

which, for our model yields, respectively 0.95 and 0.85. This 
difference could be caused by the growing importance of the 
charge transfer between 2+ and H°. According to Izotov 
et al. (2004) this effect is only noticeable when 2+ /0 + 



is larger than 8. For this model, this quotient is 7.24, and 
the relative fraction of Ne 2+ is slightly larger. On the other 
hand, the ICF(Ne 2+ ) from Mathis & Rosa (1987) is equal 
to 0.87. 



Sulphur. For this element, there are two observable stages 
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in the optical and the near IR parts of the spectrum respec- 
tively. For S+, it is usual to assume that T([Sn]) « T([Oll]), 
but in this object it is possible to measure the auroral lines 
of [Sn] at 4069,74 A. The derived [Sn] line temperature is 
consistent with that of T([On]) within the errors, which are 
large, and closer to the value predicted by the model. Thus, 
the derived and model calculated abundance of S + are sim- 
ilar. Regarding S 2+ , the line temperature agrees rather well 
with the value predicted by the model and therefore the pre- 
dicted ionic abundance is within the errors of the measured 
value. Both derived and predicted [Sill] electron tempera- 
tures agree with the value obtained from T([Om]) using the 
empirical relation derived by Hagele et al. (2006), that is: 

TQSm]) = (1.19 ± 0.08)T([Om]) - (0.32 ± 0.10) 

The ionization correction factor of sulphur accounts for 
the fraction of S 3+ in the inner parts of the nebula. The first 
assumption to calculate it was (Peimbert & Costero, 1969): 

s+ + s 2+ = o^_ 
s o 

but several authors have pointed that this aproximation 
could overstimate the fraction of S +3 in high excitation- 
nebula. Thus, Barker (1980), proposed a new ICF for S + 
and S 2+ : 



ICF(S + + S 2+ ) 



91 
O 



-l/a 



whose fit for the models by Stasihska (1978) gives a value 
of a = 3. Later models (Stasihska, 1990) have produced a 
lower value of a = 2. 

On the other side mid-IR observations for this object ex- 
ist (Nollenberg et al. 2002) and allow to measure directly the 
S 3+ abundance and calculate the ICF(S + +S 2+ ). The value 
found by Perez-Montero et al. (2006) is 1.97±0.62, which 
corresponds to Barker's formula for a value of a = 1.66. In 
that same work, a new value of the S abundance is pro- 
posed assuming a constant ratio of S 3+ /S 2+ , comparing the 
S 2+ abundances in the near- IR (9069.9532 A) and the mid- 
IR (18.71 (i). For our model, the predicted ICF(S++S 2+ ) 
is 1.31, slightly higher than the value obtained in PMD03 
(1.22) but within the errors of the new calculated value by 
Perez-Montero et al. (2006) which is 1.56±0.76. The ICF 
predicted by our model corresponds to a value of a — 2.8, 
slightly higher than the a=2.5 found in Perez-Montero et 
al. (2006). The value obtained frm Mathis & Rosa (1991) 
for this ICF is 2.28. 

Argon. Only the emission lines of [Arm] at 7137 A and of 
[Anv] at 4713,40 A are detected, thus allowing to calculate 
the ionic abundances of Ar 2+ and Ar 3+ respectively. The 
assumptions T([Arlll]) w T([Sm]) and T([Anv]) w T([Om]) 
are well justified in both cases. 

Since Ar 3+ is not usually detected in spectra, the 
ionization correction factor for argon refers to Ar 2+ . The 
ICF(Ar 2+ ) expression in Izotov & Thuan (1994) from the 
models of Stasihska (1990) yields a value of 2.51, while in 
our model a value of 1.32 is found. Using the coefficients 
from Mathis & Rosa (1991), the ICF is found to be 21.6. 
Martin-Hernandez et al. (2002) have propose the relation: 
N 2+ /N+ w Ar 2+ /Ar+. The result of our model for the first 
quotient is 7.21, while for the second one it is 41.98, what is 



fur away from the expected result. 

Iron. The only visible species in the optical spectrum of 
this galaxy is Fe 2+ through the 4658 A line. For this stage 
of iron ITL97 propose T([Fem]) « T([Oll]). In our model the 
agreement between both temperatures is good, although the 
line temperature of Fe 2+ is slightly higher. 

The ICF proposed by ITL97 from Stasihska's (1990) 
models is: 



Fe 
Fe 2 + 



1.25 ■ 



O 

o+ 



what results in a ICF(Fe 2+ ) = 10.30, slightly higher than the 
value predicted by our model, equal to 7.59, and the value 
obtained using the ICF scheme proposed by Rodriguez & 
Rubin (2004): 



ICF(Fe 2+ ) = 



o 2 + 



o 

o+ 



which gives in our model a value of 6.90. 



5 SUMMARY AND CONCLUSIONS 

We present a self-consistent model for the brightest knot 
of the Blue Compact Dwarf galaxy Markarian 209 calcu- 
lated using the photoionization code Cloudy 96. In order to 
fit the more relevant emission lines emitted by the ionized 
gas, we have introduced as input parameters some of the 
observed and deduced properties of this object, including 
radius, metallicity, WR features, the Ha flux and the rela- 
tive fluxes of the more relevant emission lines observed both 
in the optical and near-IR. The addition to the analysis of 
the far-red emission lines, including [Oll]A7319,7330 A and 
[Sm] 9069, 9532 A, allows to derive with more confidence 
the ionic abundances of oxygen and sulphur, as well as the 
ionization parameter and the equivalent effective tempera- 
ture, by means of the [Sll]/[Sm] ratio (Diaz et al., 1991) and 
the r) parameter (Vflchez & Pagel, 1988) respectively. 

Our model uses as input the spectral energy distri- 
bution produced by newly formed stellar clusters which 
has been calculated using evolutionary synthesis techniques 
(Leitherer et al. 1999) and the models atmospheres from 
Pauldrach et al. (2001) and Hillier & Miller (1998). A metal- 
licity of of 1/20 of solar (Z=0.004), as corresponding to the 
O/H abundance derived for this object, has been adopted. 
The ionization front is located at a distance of 30 pc from 
the ionizing cluster and the thickness of the ionized region is 
approximately 18 pc. The total radius is close to that mea- 
sured on Ha images. The model assumes a constant density 
of 190 cm -3 , as obtained from nebular diagnostics. 

The addition of a certain amout of dust was required 
in order to reproduce correctly the measured electron tem- 
peratures. This dust implies a gas/dust ratio of 6.13xl0~ 3 , 
identical to that in the Milky Way, and a fraction of ab- 
sorbed photons of 0.58. The dust is included in the Cloudy 
model using the command PGRAINS which provides a more 
reallistic treatment of the physics of the dust in the nebula. 
The model predicted reddening allows to correct the mea- 
sured B-R index in the region of the burst to a value in 
agreement with the predicted colour index for the stellar 
ionizing populations used in the model. Nevertheless, the 



© 2006 RAS, MNRAS 000, 1-?? 



12 E. Perez- Montero & A.I. Diaz 



logarithmic extinction at H/3 predicted by the model for the 
included amount of dust, C(H/3) = 0.19, is not compatible 
with the value deduced from the Balmer decrement (C(H/3) 
= 0.06). This could be caused by an irregular distribution 
of the dust into the gas. 

We have used the flux of Ha, the equivalent width of 
H/3 and the fluxes and equivalent widths of the WR features 
in order to constrain the age and properties of the ionizing 
cluster. Agreement between measured and predicted quan- 
tities is found for a composite population of two clusters, 
one very young, with an arbitrary age of 1.0 Myr, and an- 
other one in the WR phase with either 3.2 or 4.2 Myr. A 
marginal agreement is also found for an extended burst of 
star formation with 3.6 Myr. This scenario of continuous 
star formation is compatible with the results found by Ter- 
levich et al. (2004) from the distribution of EW(H/3) in HII 
galaxies. The masses deduced for the clusters lay in a range 
where the effect of stochastic effects in the number of mas- 
sive stars can cause variations in the predicted quantities 
but not affecting the main conclusions of our work. 

We have compared measured and model predicted elec- 
tron temperatures. The presence of the dust allows to repro- 
duce the auroral lines of [Oiii] and [Sin] and therefore the 
electron temperatures for these species, without appealing 
to any other mechanism of heating in the high-excitation 
zones of the ionized gas. Regarding the electron tempera- 
tures in the low-excitation zone, the temperature of [Sn] is 
well reproduced, while for [On], is slightly overestimated by 
the model, although it should be bore in mind that its deter- 
mination carries large observational errors. This could affect 
the derived + ionic abundance although it should not affect 
the total derived abundance of oxygen, since most of it ap- 
pears as 2+ . The electron temperatures for the rest of the 
ionic species follow the usual assumptions about the inner 
structure of the nebula and the determination of their chem- 
ical abundances are predicted to be within the errors of the 
measured quantities in most cases. The model predicts no 
fluctuations of temperature in agreement with results found 
in Blue Compact Galaxies from measurementsof the Balmer 
jump (Guseva et al., 2006; Hagele et al., 2006). 

Using the ionization correction factors predicted by the 
models for the various species, we can calculate their total 
chemical abundances and therefore we can refine the deriva- 
tion of some relevant abundance ratios. For N/O, we obtain 
-1.33, in agreement with the values expected for a certain 
amount of nitrogen produced as secondary. For S/O, the 
model predicts a value of -1.69, consistent with the constant 
value found for the sample of HII galaxies with measure- 
ments of the [Sin] lines at 9069,9532 A and which is slightly 
lower than the value measured for the Sun (Perez-Montero 
et al., 2006). This is also the case for Ne/O and Ar/O for 
which values of -0.8 and -2.34, respectively, are found. 
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